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Abstract — The Frobenius method can be used to compute 
solutions of ordinary linear differential equations by generalized 
power series. Each series converges in a circle which at least 
extends to the nearest singular point; hence exponentially fast 
inside the circle. This makes this method well suited for very- 
high-precision solutions of such equations. It is useful for this 
purpose to have prior knowledge of the behaviour of the series. 
We show that the magnitude of its coefficients can be apriori 
predicted to surprisingly high accuracy, employing a Legendre 
transformation of the WKB approximated solutions of the 
equation. 

Index Terms — Second order ODEs, Regular singular points, 
Frobenius method, Legendre transformation, WKB approxima- 
tion. 



I. Introduction 

A LARGE set of problems from technology and science 
involves the study of linear second order ordinary 
differential equations. Although the original problem is more 
likely to be a partial differential equation involving the 
Laplace operator in two or higher dimensions, it is often 
possible to reduce it to a set of one-dimensional problems 
through separation of variables. 

For example, in three dimensions the wave, heat, or 
Schrodinger equation in zero or constant potential can be 
separated in ellipsoidal coordinates (6>6)6)> related to 
cartesian coordinates by 
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where £1 > a > £2 > b > 6 > 0, plus 10 degenerate 
forms of these coordinates [ 1 1. The separated equations have 
five regular singular points, at ±a, ±6, and 00 Q. This 
means that the equations can be formulated as second order 
equations with polynomial coefficients. The majority of the 
special functions of mathematical physics, as f.i. discussed 
by Whittaker and Watson can be described as solutions 
to such equations. 

The Schrodinger equation remain separable if we add a 
potential of the form 

(6 2 - 6X6) + (6 2 - 6X6) + (6 2 - 6X(6) 
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The degenerate forms often lead to situations where two 
or more regular singular points merge to irregular singular 
points (confluent singularities). 

It is sometimes useful to evaluate solutions to much 
higher accuracy than possible with standard double-precision 
methods. Recently we have developed and used code for 
solving a large class of ordinary Frobenius type equations 
to almost arbitrary high precision, in a number of algebraic 
operations N which grows asymptotically linearly with the 
desired precision P, i.e. 

N rj iV + const P. (3) 

It was f.i. used in H to find the lowest eigenvalue of 



ip" (x) + x ijj( x ) = sip{x) 



(4) 



to an accuracy of one million decimals digits, and its eigen- 
value number 50 000 to 50 000 decimal digits. It was further 
demonstrated in [5| that it is possible to compute normaliza- 
tion integrals of the resulting wavefunctions to comparable 
precision. The method also extends to computation of many 
other types of amplitudes. 

In reference @ we have published and tested code for 
solving equations of the class 
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I.e., equations with polynomial coefficients (after multipli- 
cation by z 2 ), and at most one regular singular point in the 
finite plane. 

When implementing the Frobenius method [7 1 numerically 
the solution is represented by a convergent series 
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where the coefficients a m is generated recursively in parallel 
with the accumulation of the sum of the series. The method is 
straightforward to extend to the wider class of equations 
but this requires consideration of several special cases. 

The individual terms in |6]l may grow very big, leading to 
huge cancellations and large roundoff errors. It is therefore 
useful to have some prior knowledge of the magnitude of 
the a m 's before a high-precision evaluation — to set the 
computational precision required for a desired accuracy of 
the final result, and to estimate the time required to complete 
the computation. There is also the question whether one 
should evaluate ip(z) directly at some far away point z, 
or if it is better to make one or more steps of analytic 
continuation. I.e., evaluate ip(zi) and tp'(zi) at one or more 
intermediate points Zi. Analytic continuation of functions 



which satisfy a second order differential equation is rather 
simple to implement, since the function is fully specified 
by just two complex numbers ip(zi) and tp'(zi), plus the 
differential equation. 

We have found that \a m \ can be estimated surprisingly 
accurate from a WKB approximation of the solution, fol- 
lowed by a Legendre transform with additional corrections. 
For the general class of equations Q, or its extension to 
the WKB integrals and the Legendre transform must be 
done by numerically, but for this one can employ standard 
double-precision methods. The rest of this paper is organized 
as follows: 

In section [IT] we give a systematic presentation of all 
explicit formulas for a Frobenius series solutions of equa- 
tion (|9) around ordinary and regular singular points, also 
considering the special cases where logarithmic terms will 
or may occur. 



In section III we first give a brief motivation of the method 
of Legendre transformations, from the point of view of 
evaluating partition function integrals of statistical physics. 
We next show how the method may be improved by a 
"finite size correction", and demonstrate how it works on an 
example with known result. Although the correction method 
was motivated by our desire to estimate the coefficients of (|6| 
more accurately, it should also be applicable to the statistical 
mechanics of small systems (i.e., systems which must be 
considered before the thermodynamic limit). 

In section IV we view the sum |6]) as a partition function, 
and apply the method of the previous section to find a relation 
between the magnitudes \a m \ and IV'(-z)!- This can be used 
to estimate \a m \ provided we know 

In section [V] we discuss how one may use the WKB 
approximation to find a sufficiently good approximation to 
ip(z), thereby completing the set of tools required for our 
estimates. 

In section |VI] we demonstrate how the method works on 
a set of examples where much of the calculations may be 
performed analytically. 

A first account of this work has been presented in J8). 

II. The Frobenius method for second order ODEs 

The Frobenius method for solving homogeneous linear 
ordinary differential equations is treated in all books on 
differential equations. However, it is difficult to find general 
expressions which are sufficiently explicit for implementation 
as numerical code, and which also cover all special cases. 
We present such expressions in this section. 

We consider the second order differential operator 



C=p{z) d^ +q{z) dz +r(z) ' 



where 
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are (short) polynomials in z. Solutions to the equation 

£ip(z)=0 (9) 



in the vicinity of z = 0, depending on the behaviour of 
p(z), can often be found by series expansion. For explicit 
implementation of solution algorithms we must consider 
several cases. 

A. Ordinary points 

When po ^ the point z = is an ordinary point for 
equation ([9). The solution can be expanded in a Taylor series, 

iP(z)=J2 a ™z m . (10) 

m>0 

Insertion into equation |9]l gives 

^^P fc (m + 2)a m+2 _ fe z m = 0, (11) 



m>0 fc>0 



with 



Pk{p) = {p-k)[{fi-l- k)p k + g fe _i] + r fc _ 2 - (12) 

Here, and in the following, all coefficients with negative 
indices should be interpreted as zero, 

a- n = P-n = q-n = T- n = 0, fom=l,2 

The requirement that the coefficient of each power z m in 
( fTT| > must vanish leads to the recursion formula 
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The coefficients a and a\ can be chosen freely, or according 
to the initial conditions. 



B. Regular singular points 

When (a) po ~ with p\ ^ (and at least one of go or 
tq is nonzero), or when (b) po = p\ = qo = with p2 ^ 
(and tq is nonzero), equation |9} has a regular singular point 
at z = 0. A series solution can be found by the Frobenius 
method. I.e., one writes the tentative solution as a generalized 
power series ([6|. 

One solvability condition is that v must satisfy a second 
order algebraic equation. Hence, counting possible degen- 
eracies, there will be two solutions, v\,2 (choosing Rei/i < 
Re^). We must further consider the cases (i) v\ = v^, (ii) 
v-i = v\ + £ with £ a positive integer, and (iii) everything 
else. To avoid explicit implementation of too many special 
cases we can instead write 



ip(z) = z U2 ip{z), 



(14) 



and solve the resulting equation for ip(z). The explicit 
transformation of the polynomial coefficients is 

p(z) = z~ 1 p(z), 

q(z) = 2vz- 2 p(z) + z- 1 q(z), (15) 
r{z) = v{y — l)z^p(z) + i>z~ 2 q(z) + z~ 1 r(z), 

where p(z) and q(z) remain polynomials because po = p\ = 
<7o = 0, and f(z) is a polynomial when 



v(y - l)p 2 + vqi + r = 0. 



(16) 



I.e., when v satisfies the indicial equation. The transformed 
equation corresponds to the case (a) above, with index = 
and qo/pi > 1. The latter implies v\ < 0. Henceforth we 



will only consider this case, dropping " from the notation. 
For the computations below we need various cases of the 
formula (using Po(/i) = when p = 0), 



£ 52 ( a °< m + ai ' m log z ) z * 

m>0 

52 [-Pfcfi (m+2+u) (a ,, n +i-fe + a^mfi-fc log z) 

m>-i k>a 

+ Q m (m+2+v) ai,„ l+ i_ fc ] }z m+ ", (17) 
with Pfe(/x) given by equation ( fl2| ), and 

Q fc 0*) = (2/i-l-2fc)p fc +gfc_i. (18) 

C. Regular singular point with non-integer index difference 

First consider the everything else case. We assume a 
solution of the form j6l. Now Cij){z) — implies 



52 Pk+i{m + 2 + i/) a„ l+ i_ fc = 0, 



fc>0 



for m = — 1, 0, . . .. For m = — 1 this becomes 
Pl(l + v) = u[(y-l)p 1 +q o ]=0, 



(19) 



(20) 



with solutions and V\ = 1 — qo/pi, and iaj = 0. In this case 
2/1 is non-integer. For m > 0, and i/ = V\ and £/ = 1/2 = 



respectively, equation (19i implies the recursion relations 
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(m+l)(m+l — i/i)p; 



k>l 



In both cases the coefficient clq can be chosen freely, or 
according to initial conditions. The last recursion in (21 1 is 



always working, since by arrangement v\ < 0. The first one 
is working in this case since 1/1 is non-integer by assumption. 



D. Regular singular point with degenerate indices 

For the degenerate case of v\ = v-z = we still have a 



solution corresponding to the last recursion in equation (21 
For a linear independent solution we make the ansatz 



%/}(z) = 2J ( a 0,m + ai. m log z) Z T ' 
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Now Cijj(z) = implies 

52 P fc+ i(m+2) ai )JTtf i_ fc = 0, 



fc>0 



(23) 



52 Pfe+i(™+ 2 ) ao,m+i-fc + Qfefi(m+2) a ljm+1 _ fc = 0. 
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for m = — 1,0, For m 

Pi(l) = 



-1 this becomes 
?i(l) = 0, 



(24) 



which holds since go = Pi when v\ = 0. For m > 



equation ( 23 1 implies the recursion relations 
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The coefficients ao.o and a^o can be chosen freely, or 
according to initial conditions. 



E. Regular singular point with integer index difference 

For the case v\ = —I (a negative integer) and 1/2 = we 
still have a (1/2) solution corresponding to the last recursion 



in equation (21 1. For a linear independent solution we make 
the ansatz 



1p(z) = 52 ( a 0,m + ai.ro l°g z ) • 
m>0 

Now Ctp(z) = implies 

52 Pfcfi( TO + 2 -^) ai,ro+i-fc = 0, 

fc>0 



(26) 



(27) 

52 ^W m + 2- '0 a o,nH-i-k+Qk+i(m+2-£) ai jm+ i_fc = 0. 

fc>0 

for m = —1,0,.... For m = — 1 this can be solved by 
chosing a 0j o arbitrary, since Pi(l— I) = when V\ — —I and 
q = (l+£)pi, and a 1)0 = since Qi(l-£) = -^pi + 1 0. We 
may next use the first recursion in equation pi) to compute 
ao,i+m until m+1 = £, where it breaks down. For m+1 = I 
equation (|27| becomes 
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(29) 



52 ^+i( 1 ) a o/-fe + Qi(l)av = 0, 
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with solution (note that Qi(l) = £pi) 



«i, 



^P fc+1 (l)a , 
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We may choose ao ^ freely; this correponds to an addition 



of the i/ 2 solution. For m + 1 > I equation (27 1 implies the 
recursion relations 
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This completes the description of the Frobenius series so- 
lution of Cip(z) = around ordinary and regular singular 
points. 



III. Legendre transform with finite size 

CORRECTIONS 

Although the Legendre transform is a well-defined math- 
ematical procedure by itself |9|, 1 1()|. ifTTl . it can f.i. be 
motivated as a maximum integrand approximation to integral 
relations between partition functions in equilibrium statistical 
physics. As such it can be modified to take into account 
corrections due to (gaussian) integral contributions around 
the maximum. 



A. Maximum integrand approximation 
Consider the integral relation 



dye 



~U(y)+p-y 



(31) 



where U(y) is assumed to be a smooth function of its 
argument. We first approximate the integral by the maximum 
value of the integrand. This occur at some point x where 

-V x U{x)+p = 0, 

and the approximate value of F(p), which we shall denote 
Fq(p), becomes p-x — U(x). These relations correspond to 
the standard Legendre transform from the pair x, U(x) to 



p=V x U(x), 
Fo{p) =p-x-U(x) 



(32) 
(33) 



where x is assumed to be eliminated by use of equation ( 32 1. 

Consider next the effect of changing p by a small amount, 
p — >• p + dp. This will change the maximum a; by a small 
amount, x — > x + Sx. This leads to a small change in Fq, 

F (p) ^F (p + Sp) = F (p) + Sp ■ V p F (p) 
= (p + Sp) ■ (x + Sx) — U(x + Sx) 
= p ■ x — U(x) + Sx ■ [p — VajE/(a:)] + Sp ■ x, 



when neglecting second order corrections. Eliminating terms 
by use of (32i and (33 i gives Sp ■ V p Fq(p) = Sp ■ x, or 



since the direction of Sp is arbitrary, 

x = V p F (p), 
U{x) = p ■ x - F (p). 



(34) 
(35) 



This is the inverse Legendre transform. Equation ( 35 i is just 



a trivial rewriting of ( 33 i 



It follows from equation (32i that the Jacobi matrix 
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U{x), 



and from equation (34i that the Jacobi matrix 
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■F (p). 



Since Jacobi matrices of inverse transformations are matrix 
inverses, we have the relation 
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B. One loop fluctuation correction 

We next improve the evaluation of (31 1 by expanding U(y) 
around its minimum. We write y = x+z and assume that the 
main contribution to the integral comes from a small range 
of z-values. Hence 

U(x + je) « U{x) + z ■ V x U{x) + izV'[/ y -(a:), 

which gives 

e F(p) =e p- x -U( x ) J ^g-^VtTij)*) 

==e P-»-ir(«) det{U l . J {x)/2Tiy 1/2 

I.e., denoting the corrected expression Fi(p), 
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Fi{p) = px 
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-logdet(U^(a;)/27r) 



F (p)- -logdet(tf -(a!)/27r) 



(38) 



This corresponds to a standard one-loop correction to the 
partition function e F ' p ) of statistical physics, where one 
may proceed with an ordinary inverse Legendre transform 
to define a (fluctuation corrected) effective potential 



x = V p Fiip), 
r(a;) = p ■ x — Flip). 



(39) 
(40) 



However, here we want to compute the original potential 
Uix) from the information provided by Flip). The rela- 
tion d36ll now reads 



Uij{x) 



d 2 F Q (p) 



(41) 



dp" 1 dpi 

but we have no direct access to the quantity Fq(p). Equa- 
tion ( [38) instead becomes a second order differential equation 
relating Fq(p); hence an exact inversion seems difficult in 
general. However, if the one-loop correction is assumed to 
be small compared to Fo(p), we may approximate 

'd 2 F ip)\ (d 2 Fiip)\ _ 



and 
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logdet(2 7 rF 1 . iJ (p)). 



(42) 



We finally apply the inverse Legendre transform ([34]), ( |35] l 
of the pair p, -Fb(p) to find the pair x, Uix). 

C. Simple application 

Consider the "partition function" 

e'W = 2~ N (1 + e p ) N . (43) 

This actually is the generating function for the probability 
of x heads in a sequence of N independently flipped coins, 
corresponding to the probability 



P x = e- u ^ = 2 



—N 



N 



(44) 



but assume we don't know that. We instead use the method 
above, with 

Hp) = Anog(l+e p )-An g2, (45) 
f'ip) = Ne p il+e p )- 2 . (46) 



Approximations to the binomial distribution for N = 10 



— Legendre approximated distribution 
o o o Stirling approximated distribution 
• • • Exact binomial distribution 




Approximations to the binomial distribution for N = 50 




Fig. 1. Comparison of the continuous distribution obtained by a Legendre 
transform (with finite size corrections) of the generating function, cf. equa- 
tion (48) , with the exact binomial distribution {44} for N = 10. For a fair 
com parison one should perhaps com pare the exact probability P x given by 
{44| with the average value of {48| over the interval (x — + The 
result of making a Stirling approximation of the binomial coefficient is also 
shown. This approximation gives P x = at the endpoints x = and 
x = N. 



From 



dp 

u(x) = px — 



/(p)--log27r/» 
/(p)-ilo g 27r/ , '(p) 



we find 



p = log£ - log(l - £), 



where we have written x = (TV + 1)£ — ^ to simplify 
expressions, and 

u(x) = (N + 1) K log£ + (1 - log(l - 0} 

+ |log27riV + iVlog2. (47) 

This means that the discrete probability P x of equation (44 1 
is approximated by a continuous distribution 

1 



-u(x) _ 2~ N 
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(48) 



which should be compared with the Stirling approximation 
to the binomial coefficient 



N 



1 



-N[ X log X +(l-x)log(l-x)] 



(49) 



where x = N\- Comparisons between the exact distribu 



tion (44 1, its Stirling approximation (49 1, and the Legendre 
transformed approximation (|4"8j ) are shown in Fig. [T] and [2] 

IV. Estimating Frobenius coefficients by 
Legendre transform 

We now use the results of the previous section for our 
main objective, to estimate the magnitude of the coefficients 
a m . Here we will not consider the presence of logarithmic 
terms. Our basic hypothesis is that the sum (|6ji for large \z\ 



Fig. 2. Comparison of the continuous distribution obtained by a Legendre 
transform (with finite size corrections) of the generating function, cf. equa- 
tion (48} , with the exact binomial distribution {44) for N = 50. The result 
of making a Stirling approximation of the binomial coefficient is also shown. 
This approximation gives P x = at the endpoints x = and x = N. 



receives its main contribution from a relatively small range 
of m-values, at least for some phase values of z. Introduce 
quantities u and s(m) so that 



Hence our hypothesis is that 

= maxV>(e" +i! : 



l s(m) + (i/+m)u 



(50) 



with the main contribution to the sum coming from a 
relatively small range of m-values around a maximum value 
m. The latter is defined so that s' (rh) + u = 0, s" (m) < 0. 
Now write m = ffi + Am, and approximate the sum ( ]50] > 
over Am by a gaussian integral. This gives 

e 5 ^ a v /-2 7 r/s"(m)e s ^ + ^ + ™) u . 

In summary, we have found the relations 

u = — s (m), 



(51) 



S{u) = s(m) — (y + m)s'(m) + — log 



2tt 



'(m) 



S {u) + - log 
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-s"(m) 



(52) 



Compared with the discussion in the previous section we see 
that this is essentially a Legendre transformation between 
s(m) and S(u), but since there are some minor differences 
we repeat the derivation of the inverse transformation. 

Consider a small change u — >• u + Su. To maintain the 
maximum condition we must also make a small change m — > 
m + 8m, with 8m = —8u/s"(m). I.e. s"(m) — —u'(m). 
This is consistent with the result of taking the m-derivative 



of equation (51 1. One further finds that Sq(u) becomes 



S (it + 8u) = S (u) + S' (u) 5u + ~S'^{u) Su 2 + 



s(m) + [y + m)u + (m + u) 8u — 



1 



2s" (m) 



8u 2 



giving the relations 

(m + u) = S' (u), 

s{m) = Sq(u) — uS' (u), 
s"( m ) = -S'oiu)- 1 . 



Numerical and estimated coefficients 



(53) 
(54) 
(55) 



Equation (55 i just says that (dm/du) = [du/ dm)' 1 . We are 
only able to compute S(u) directly, not Sq(u), However, they 
only differ by a logarithmic term, hence we will approximate 

log(-s"(m)) = - log S%(u) w -logS"(u). This gives 

So («)«£(«)-£ log (2tt £''(«)), (56) 



which can be used in equations ( |53] - [55] l when we have 
computed S(u). 

V. WKB APPROXIMATION 

It remains to find S(u). Here we will use the leading order 
WKB approximation to find a sufficiently accurate estimate. 
When z = is an ordinary point, i.e. when i/_ = 0, v + = 1, 
the leading order WKB solution to (g]> is Il2l. lull lfl4l 

VQo/Q(z)exp(lJ (57) 

where Q 2 (z) — Yl n =i v n zn ~ 1 < an ^ Qo — Q(0). This repre- 
sents a superposition of the solutions ip±(z). The difference 
between the v + and i/_ solutions is at worst comparable to 
accuracy of our approximation; hence we will not distinguish 
between them. 

When z = is a regular singular point we use the 
Langer connected WKB approximation to obtain leading 
order solutions in the form 



z u± VQo/Q(z) x 
1 f z At 

T 



exp ± 



(58) 



Here Q 2 (z) 



1 »2 



E 



_n+l 



and Qo 



Q(0). In equation (58 i we distinguish between the v+- and 
i/_ -solutions, because the difference v+— v- may in principle 
be large. 

The WKB integrals must in general be done numeri- 
cally, sometimes along curves in the complex plane. This 
requires careful attention to branch cuts. One must also 
take into account that the behaviour of the WKB solution 
may change from one exponential behaviour to another, due 
to the existence of Stokes lines [15|. We have observed 
that this sometimes can be explained as contributions from 
topologically different integrations paths. In this paper we 
will only give some examples where most of the calculations 
can be done analytically. 

VI. Examples 
A. Anharmonic oscillators 
Consider the equation 

for real c so that c 2 > 0. For large y the typical solution 
behaves like 

~ e ^ 3 + c \ (60) 
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Fig. 3. Comparison of nu merical coef ficie nts a m (points) with estimates 
(full-drawn lines) based on j62||63| and {74} . The estimates of log \a m \ are 
accurate up to corrections which depend logarithmically on m. 



neglecting the slowly varying prefactor. For a given value of 
\y\ this is maximum along the positive real axis. Hence, with 
x = y 2 = e 11 , we find as a leading approximation 

S(u) = | (e5" + 3c 2 e2" 

In this case the Frobenius series can be written 



m=Q 



m=0 



(61) 



(62) 
(63) 



with v = 0, 1. Ignoring the log(S"(u))-term in ( [53] [54} [56 
we find 

log (\a m \) = (| - \u) el" + c 2 (1 - \u) e* u . 

For c = an explicit representation is 

2 

log |a m | = -m(l-log2m). (64) 

This is plotted as the lower curve in figure [3] It fits 
satisfactory with the high-precision coefficients generated 
numerically, but there remains a correction which depends 
logarithmically on m. For nonzero c the parametric represen- 
tation provides equally good results, as shown by the upper 
curve in figure [3] 

The conclusion of this example is that for a fixed (large) 
x we expect the largest term of the power series to be 



max|A m (a;)| - ea ( 



(65) 



neglecting a slowly varying prefactor. Further, the maximum 
should occur at 



i (V/ 2 + cV/ 2 ) 



(66) 



Finally, estimates like equation (64i for the coefficients a m 
may be used to predict how many terms A4 we must sum to 
evaluate *p(x) to a given precision P, based on the stopping 
criterium 

(67) 



\a M \x M < IQ- P . 



As can be seen in figure |4] the agreement with the actual 
number of terms used by our evaluation routine is good, 
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Fig. 4. This figure compares the a priori prediction, based on equation (64) , 
of the number of terms Ai which must be summed in order to evaluate 
^(y) for c = to a desired precision P, with the actual number of terms 
computed by the numerical routine 1 6 1 . 



Fig. 5. This figure shows the ratio between the computed coefficients 

a m and the crude pr edic tion 464i (labelled la^ I) and the logarithmically 

I I — (l) 
corrected prediction |70l (labelled \arn |). For easy comparison we have 

in both cases adjusted an overall constant such that the ratio is unity for 

m = 3. 



in particular for high precision P. But keep in mind that a 
logarithmic scale makes it easier for a comparison to look 
good. 

Next consider the logarithmic corrections. Including the 



prefactor of equation ( 60 » changes S(u) by an amount 
AS-(u) = -ilog (e" + c 2 ). 

Including the log(S"'(u))-term in the relation between S(u) 
and Sq(u) changes Sq by an additional amount 



AS (u) = -ilogQe§" + i C 2 e^ 



(69) 



For c = this changes the relation ( 64 1 to 



log M = I (2m + 5/2) (1 - log (2m + 5/2)) . (70) 



For \a m \ this essentially corresponds to a factor m 



-5/6 



B. Double well oscillators 

The same procedure also work for the equation 

-j^*(s/) + (2/ 2 - c2 ) 2 *(y) = °> < 71 > 

which however is a little more challenging since the maxi- 
mum value of ^(ye'^l sometimes occur for ip =^ 0, i.e. for 
complex arguments. 

For large y the typical solution behaves like 



• es ! 



(72) 



neglecting the slowly varying prefactor. Equation f7T| ) can be 
transformed to the form (pb by introducing x — y\^{y) = 
ip{x). Hence, with x = y = e" 



S(u) 



max 5 

v 6 



|Re ( e i(" +i ^ 



3c 2 e |(«+V) 



1/2 

The maximum occurs for cos \y> = (l + c 2 e _ ") 
when e 11 > ^c 2 , and for cos \tp = —1 otherwise. This gives 

„2 C 



c 2 e «/2 _ l e 3u/2 for e u < %c . 



±C 2 

v ' 1 i(e n + c 2 ) 3 / 2 for e u > \c 2 



This implies that 



ie u / 2 (c 2 - e 11 ) fore u <|c 2 , 
|e u (e u + c 2 ) 1/2 for e" > |c 2 , 



(74) 



(68) log(|am|) 



(1- \u) c 2 e"/ 2 -(!-!„) e 3 "/ 2 for e" < |c 2 , 
[(|-|«)e" + |c 2 ] (e u + c 2 ) 1/2 for e" > |c 2 . 



This representation compares fairly well with the numeri- 
cally generated coefficients, as shown by the middle curve 
in figure [3] However, in this case the coefficients a m have 
a local oscillating behaviour. The representation ( 74 1 should 
be interpreted as the local amplitude of this oscillation. 

The conclusion of this example is that we expect the 
largest term of the power series to be term of the series 
to be 

maxK^I-e^ 2 ) 3 ' 2 , (75) 

m 

neglecting the slowly varying prefactor. Further, the maxi- 
mum should occur at 



\x (x + c 2 ) 1/2 ~ iz 3 / 2 + icV/ 2 . 



(76) 



VII. Conclusion 



As illustrated in this contribution the coefficients of Frobe- 
nius series can be predicted to surprisingly high accuracy 
by use of Legendre transformations and lowest order WKB 
approximations. We have also tested the validity of the 
method on many other cases. 
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